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On the Puzzle of Odd- Frequency Superconductivity 
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Since the first theoretical proposal by Berezinskii, an odd-frequency superconductivity has encountered 
the fundamental problems on its thermodynamic stability and rigidity of a homogenous state accompanied 
by unphysical Meissner effect. Recently, Solenov et al. [Phys. Rev. B 79 (2009) 132502.] have asserted 
that the path-integral formulation gets rid of the difficulties leading to a stable homogenous phase with an 
ordinary Meissner effect. Here, we show that it is crucial to choose the appropriate saddle-point solution 
that minimizes the effective free energy, which was assumed implicitly in the work by Solenov and co- 
workers. We exhibit the path-integral framework for the odd-frequency superconductivity with general 
type of pairings, including an argument on the retarded functions via the analytic continuation to the real 
axis. 

KEYWORDS: odd-frequency superconductivity, Stratonovich-Hubbard transformation, saddle-point approx- 
imation, Meissner kernel 



1. Introduction 

The field of anisotropic superconductivity has developed 
extensively in condensed matter physics. In the anisotropic 
superconductors, a strong Coulomb repulsion with relatively 
local character tends to suppress an on-site amplitude of the 
Cooper-pair wavefunction, and it favors an off-site pairing in- 
stead. It is also possible to form an off-time or equivalently 
an odd-frequency pairing to avoid the strong repulsion. It was 
first proposed by Berezinskii in the context of the superfluid 
3 He, [ ' and was revisited later for the mechanism of the cuprate 
superconductors by Balatsky and Abrahams. 2 ' After the pro- 
posal of the odd-frequency superconductivity, considerable 
amount of investigations have been carried out both theoreti- 
cally 3-10 ' and experimentally 11 ' in a wide variety of theoreti- 
cal models and realistic materials. The anomalous proximity 
effect in the superconducting junctions and superfluid 3 He has 
also been investigated by the odd-frequency pairing mecha- 
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In contrast to proceeding of such stimulating investigations, 
fundamental problems lying in the odd-frequency supercon- 
ducting state remain unresolved. The nuclei of the prob- 
lems are thermodynamic instability of the homogeneous odd- 
frequency superconducting state 5 ' 11 ' 18 ' and unphysical Meiss- 
ner effect with the negative Meissner kernel. 6 ' In other words, 
these issues end up with a pure imaginary penetration depth 
and a deficiency of rigidity of a homogeneous order parame- 
ter, and that the free energy of the superconducting phase is 
lower than that of the normal phase in the high-temperature 
limit. The common source of these problems is ascribed to 
the relation among the anomalous green functions, 



F + (i(i)„) - F(-ia>„Y = -F(iw„y, 



(1) 



in the case of the odd-frequency pairing (see the discussions 
below). In the case of the even-frequency pairing, the opposite 
sign appears in the most right-hand side, and no problems as 
mentioned in the above do occur. 
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Recently, Solenov et al. have demonstrated that the path- 
integral treatment of the odd-frequency superconductivity 
ends up with the different relation from (1), i.e., 

F + {ioj n ) = F(ioj n Y, (2) 

which is the same relation for the even-frequency pairing. 19 ' 
As a result, a homogeneous odd-frequency superconductiv- 
ity is fhermodynamically stable in the bulk and exhibits the 
ordinary Meissner effect. In their argument, the lack of an ap- 
propriate hamiltonian with broken U(l) gauge symmetry is 
vital in the case of the odd-frequency pairing, and only the 
path-integral approach provides reasonable results. Although 
this is an important and correct statement, the essence of the 
problem still seems to remain obscure. 

In this paper, we aim to unveil the essence of the prob- 
lems by emphasizing the choice of the saddle-point solution 
for the odd-frequency pairing. Although the correct choice 
between the gap function A(iu>„) and its conjugate counter- 
part A(ioj„) is not apparent, Solenov et al. have implicitly 
assumed A(ia>„) = A(/w„)* in their argument [eq. (10) in 
ref. 19]. In view of the delicate issue to remedy the fate of 
the odd-frequency superconducting state, we should demon- 
strate a detailed formulation for the odd-frequency pairing 
state with great care in the path-integral framework. 20 ' 21 ' In 
the next section, we give the derivation of the free energy, 
and discuss the correct choice of the saddle-point solution for 
the case of an odd-frequency equal spin pairing as the sim- 
plest example. With the correct choice of the saddle point, 
we show that the Meissner kernel becomes positive as similar 
to the even-frequency pairing. The green functions and then- 
analytic continuation to the real axis are also given. We sum- 
marize the paper in §3. The formulation for the general type 
of pairings is given in Appendix. 

2. Path-integral framework for superconductivity 

In this section, we derive the minimal formula to argue the 
cure for the odd-frequency superconductivity. In order to fo- 
cus on the essence of the issue, we consider the case of the 
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triplet pairing only with f spin pairs, i.e., the pairing interac- 
tion V (< 0) is assumed to work only among f spins (it is 
equivalent to the case of the spinless fermions). 

2.1 Free energy 

We first give the derivation of the effective free energy func- 
tional for superconductivity. Let us begin with the partition 
function in the path-integral form, 

Z = J Dt//Diffe~ (So+Sint) , 

S = J^ t (l)(3 T +£Wt(l). 

S-m = ^ JJV(l-2)p n (l,2)p n (l,2), (3) 

where we have used the abbreviations for the space-time coor- 
dinate, 1 = x\ = (ri,n) and jj = f dn f dr\ for notational 
simplicity. /3 = l/T is the inverse temperature, and the ki- 
netic energy, £ = -V /2m - p is measured from the chemical 
potential, p. The pairing interaction V(l - 2) is a real even 
function with respect to {x\ - x^). A strongly retarded inter- 
action may be important to realize the odd-frequency pair- 
ing. p n (l,2) = i/^ T (l)(A T (2) andp n (l,2) = t// T (2)t// T (l) are 
the pair-density field and its complex-conjugate counterpart. 
The anti commutation (Pauli principle) demands the relations, 
p n (l,2) = -p TT (2, 1) andp n (l,2) = -p n (2, 1). 

Using the Stratonovich-Hubbard transformation for the in- 
teraction part, 20) we have 

S aux = ~ ||[p Tt (l,2)A Tt (l,2) + A n (U)p n (l,2)], 

5a= 4IIw^)^ (1 ' 2)Att(1 ' 2) - (4) 

Note that the auxiliary fields, A T1 -(1,2) and A T1 -(1,2), are the 
complex onumber fields corresponding to the gap functions, 
and they satisfy the same symmetry relations of pn(l, 2) and 
p TT (l,2). At this point, we have introduced no approxima- 
tions. 

Now, let us find a homogenous saddle-point solution that is 
independent of the center-of-mass coordinate, R = (x\+X2)/2. 
For this purpose, we replace the path integral over A TT (1,2) 
and Aft (1,2) with trial path for the saddle point. Then, we 
define the "mean-field" free-energy functional as 



/3T M f[A, A] = - In Z MF = - In 



DifrDif/e 



-SmfGM'AA) 



(5) 

Here Smf = Sq + S& + S dux is the "mean-field" action in 
which Af|(l,2) and A T1 -(1,2) are replaced by the trial path, 
A-|-f(r) and Ajj(r), with the relative coordinate r = x\ - X2. 
The true saddle point will be determined by minimizing the 
free-energy functional with respect to the trial path, A| T (r) 
and A-f|(r). 

Let us introduce the fourier transformation, 



= -= X ukv kxi . w = -j* 2 H k y~ ,! ' 



A TT (r) = A n (k)e ikr , A TT (r) = A n(k)e~ ikr , 



(6) 



where k = (k,ia>„), tj n = (2n + l)n/fi, 2Zk = Hi Xm an d 
kx = k ■ r - co n T. Then, the action containing the fermion 
fields is expressed as 



So +S, UX = ^J] (W> M-k))M u (k) 

1 k 

with 



ffft(-k))' 



= 1 ArtW 



A n (k) 
-ioj n - %- k 



(7) 



(8) 



Integrating out the fermion fields in (5), we obtain the explicit 
expression of the "mean-field" free-energy functional (assum- 
ing = 



Tmf[A, A] 



2/3 f r 



A TT (r)A TT (r) 
V(r) 



IVln 

2/3^ 



detM TT (fc) 



detM TT (fc) = -[co 2 n +f k + A TT (*)A n (*)], 



(9) 



where we have measured the free energy from that in the nor- 
mal state. Note that since we have used no explicit relations 
for the odd-frequency pairing, the obtained free-energy func- 
tional is also valid for the even-frequency pairing as well. 

2.2 The correct choice of the saddle-point solution 

In the path integral, (4), the contours for Afj(l,2) and 
Aff(l,2) can be taken independently. 2 ^ However, in the 
saddle-point approximation, Ajj(r) and A^(r) are no longer 
independent with each other. Taking account of the fact that 
the free energy must be real, we have two possible choices of 
the extrema. Namely, we have the choices, 22) 

A TT (r) = +A TT (r)\ A TT (/fc) = +A TT (fc)*, (10) 

or 

A TT (r) = -A TT (r)\ A n (fc) = -A n (jfc)*. (11) 

In order to demonstrate the difference between these two 
choices, let us consider a gap function of the form, Aff(r) = 
A () sin(7TT/y6) with a real Ao (an odd-frequency pairing) for 
example. After some manipulation, we obtain the mean-field 
free energy near the transition temperature, T c as 

Tmf = ±a (|- - l) A 2 + bA 4 + O(A 6 ), (12) 

where a and b are positive constants. The upper or the lower 
signs correspond to the choice (10) or (11). From this con- 
sideration, we must choose the choice (10) for a saddle -point 
solution, yielding Ajj(r) = Ajj(r) = in the high-T normal 
phase, while A^(r), A^(r) + in the low-T superconducting 
phase. Note that the relation (10) is implicitly assumed in the 
work by Solenov et al. [eq. (10) in ref. 19]. 

Together with the condition (10), the true saddle-point path 
is given by minimizing the free-energy functional (5) as 



o=/^ MF 



6S 



MF 



5A TT (r) WA TT (r)/ MF 







MF \ 

^TtW/mf' 



(13) 
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where we have defined the statistical average with respect to 
the mean-field, 



(■ ' - )mf - 



•SmfGM-AA) _ 



(14) 



Equation (13) gives nothing but the gap equations, 

A TT (r) = -V(r)F n (r), A TT (r) = -V(r)F n (r), (15) 

where we have defined the anomalous green functions, 

f> T (r) = <p TT (l,2)) MF , F TT (r) = (p n (l,2)) MF . (16) 

Once we choose the correct saddle point by (10) and (15), 
it is easily confirmed with the help of the gap equations (15) 
that the associated anomalous green functions (16) satisfy the 
relation, 

F n (r) = +F n (ry, F n (k) = +F tT (*)*, (17) 
where the fourier transforms of F^(r) and F^(r) are given by 

F n (k) = <<M*Wr t (-*)>MF • F Tt» = (^t(-£¥T W) MF ' 

(18) 

We can also confirm this relation explicitly from the definition 
as was shown by Solenov and co-workers. 19 ' 

Now, let us make some comment on the relation between 
F^(k) and F^(k). For simplicity, we consider an on-site 
s-wave pairing since the frequency ia> n dependence of the 
anomalous green functions is essential here. If we were using 
the hamiltonian formalism, we could write down the spectral 
representation of the anomalous green functions as 

F n do>n) = ~zh jz +E - f \ n r t I m ) \ m K r) ' 

" run " " 



F n (i(o n ) = ~Yj 



Em 

-f3E m + e -f3E„ 



Zi idjfj ~\~ -E pi Eif) 



(n | if/^ | nij (m | i/^ | nj , 
(19) 



where Z = e~^ E ". From these expressions, we would im- 
mediately obtain the relation F^(ico n ) = F^(-ia> n y . Then, 
the property of the odd-frequency pairing would eventually 
lead to Fii(iu>„) = -F^{ia> n )*, which seems to contradict 
with (17). As was pointed out by Solenov et al., 19> however, 
there exists no appropriate mean-field hamiltonian, Hmf in 
the presence of the retardation of the gap function, although 
the above argument is based on the relation, H MF \n) = E n \n) 
and (n\ H MF = (n\ E n . Therefore, the above proof cannot work 
out due to the explicit use of the mean-field hamiltonian. Note 
also that if we chose the maximum solution of the saddle- 
point (1 1), we would obtain the incorrect relation, F^(uo„) = 
-Ftf(i(o n y. It is the relation that has been widely believed to 
hold in the case of the odd-frequency pairing. It should be 
also emphasized that even for the even-frequency pairing the 
ordinary proof for F^(ia) n ) = F^(-icj n )* = F^(io) n )* with 
the explicit use of Hmf should not work in the presence of 
the retardation. Unfortunately, we did not realize the potential 
problem on this aspect since in the case of the even-frequency 
pairing the relation in question is continuously connected to 
the one with no retardation effect where the proof with H MF 
should work out. It only becomes evident when we have con- 
sidered the odd-frequency pairing. 

The discussions so far made in literatures (except 



ref. 19) have begun with the incorrect relation, F^(ioj n ) = 
-Ftf(ia> n y, and the relation between A^(k) and A^(k) are 
determined afterward so as to be consistent with the incor- 
rect relation. This is the source of various drawbacks of the 
odd-frequency pairing. However, since the anomalous green 
function should be defined with respect to the given mean- 
field gap function with broken U(l) gauge symmetry, we must 
first set up the relation between A^(k) and A^(k) by the min- 
imum condition of the free energy as (10), and the anomalous 
green functions are determined afterward. With this proce- 
dure, we always obtain A^(k) = A^(k)* and F^(k) = F^(k)* 
both for the even-frequency and the odd-frequency pairing. In 
summary, there are no formal differences in various formula 
between the even- and the odd-frequency pairings except the 
symmetry in the frequency dependences. 

2.3 Meissner kernel 

Let us consider the magnetic-field response. For simplicity, 
we consider the London limit where a spatial dependence in 
A can be neglected. With the London gauge V • A = 0, the 
effect of the magnetic field is taken account by replacing the 
gauge-invariant minimal coupling as V — > V + (ie/c)A in £ 
(e > 0). Namely, we replace the matrix M^(k) with 



Cm n ^ Cm e (2ck ■ A + eA 2 

M n (k;A) = M n (k) + —^ 









2ck ■ A - 



eA 1 )' 
(20) 



Then, the expansion of the free energy up to 0(A ) yields 

2 W 



Tmf(A) 



2 mc 1 4^ 



(21) 



>j 



where the Meissner kernel has been introduced by 



K u (0) 



-J] 



(G TT (fe) 2 + F n (*)F n (*)*) + 6ijGf[(k) 
m 1 ' 



The current density is then obtained as 
dTuF(A) 



Ji 



dAi 



= — yK u (Q)Aj. 

mc /_j 



(22) 



(23) 



Note that we have used the relations (10) and (17). Here, the 
explicit form of the green functions is given by 

-(ioj n + &) 



G u (k) = 
F n (k) = 



<jg+$ + |A n (*)p' 
A TT (fc) 



coi+e k +\A n (k)\ 2 ' 

which are obtained by the matrix inverse as 



-G^(k) Fff(k) 
F n (ky G n (-k) 



M^(k) 



-\ 



(24) 



(25) 



Since the above expressions are exactly the same as the con- 
ventional ones 23) even in the case of the odd-frequency pair- 
ing, the Meissner kernel (and the superfluid density) is posi- 
tive as the consequence of the relation (10). 

2.4 Analytic continuation to the real axis 

Owing to the general analytic property of the green func- 
tions, the gap function can be expressed in the form (omit the 
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index k for simplicity), 

j w A' T ' T (w) 



A TT (z) 



J-co 7!" 



10-Z 



, A TT (z) 



J-oo 



-, (26) 



where A^(w) = Im A|-f(w + /0) and A^(w) = Im A^(w + iO). 
By the relation Afr(-Ki>„) = <pA^(iaj„) (0 = +1: the even- 
frequency pairing, (p = - 1 : the odd-frequency pairing), we 
have 

A" T (-w) = -<£A" T (w), A" T (-w) = -0A" T (w). (27) 

From this relation, A^(iu>„) and A|-f(/w„) is real (pure imag- 
inary) for the even-frequency (the odd-frequency) pairing. It 
is also shown that the retarded functions satisfy the relation, 



Afi(—cj) — <pA^(a>)*, Ai"f(— a>) — 0Aff(a>)* 



(28) 



Thus, the terms "even-frequency" and "odd-frequency" rep- 
resent the symmetry property of the Matsubara and the real 
part of the retarded gap functions. 

Moreover, the minimum-energy condition (10) demands 



A^(oj) — (pA'^(oj), 
and then the retarded functions satisfy 
A|f(o>) = <pAfi(to). 



(29) 



(30) 



Therefore, the analytic continuation from the upper half plane 
to the real axis is carried out by the replacement, 

Aff(zw n ) — > Aff (w), A^(ico„) = A^(i(o„Y — > (/>A^(oj). 

(31) 

Then, the retarded green functions are given by 



F n (k, oj) 



D 1 - % 2 k - 0A T f(£, to) 2 + /Osgn(w) ' 

- A n (k, u) 
to 2 - %\ - <pA^(k, cS} 2 + /Osgn(w) 



(32) 



Here the global phase of the gap function due to broken U(l) 
gauge symmetry does not appear in the denominator of the 
green functions. Note that the retarded anomalous green func- 
tions satisfy the similar relations of the gap functions, and 
<pA(k, oj = 0) 2 > 0. 

3. Summary 

We have derived the free energy of the superconducting 
state as a saddle-point solution in the path-integral frame- 
work. The emphasis has been placed on the correct choice of 
the minimum free-energy solution that provides us with the 
ordinary relation between the gap function and its conjugate 
counterpart, 

Al a {k) = A afs (k) = A al3 (ky. 

The associated anomalous green functions determined in a 
consistent way through the gap equations satisfy the similar 
relation, 

^(t) = ^(t) = W. 

irrespective of the symmetry of the gap functions in the fre- 
quency domain. A naive derivation for this relation using the 
mean-field hamiltonian should not be valid in the presence of 
the retardation in the gap functions. 

The gap function in Matsubara representation is pure imag- 



inary for the odd-frequency pairing. The analytic continuation 
from the upper half plane to the real axis is carried out by the 
replacement, 



A aP {k, ico n ) -> A a/3 (k, u), 
A' (£, ioj n ) = A a/j (k, ia>„y 



-A afi (k,cj), (33) 



with A a p(k, -co) = -A a p{k, cj)* for the odd-frequency pairing. 

From now on, the odd-frequency pairing has no apparent 
deficiency, and hence a proper strongly retarded interaction 
can mediate the odd-frequency pairing in real materials. 
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Appendix: In the case of general type of pairings 

In this appendix, we repeat the discussions in the case of 
general type of pairings. The formula obtained here are valid 
both for the even- and the odd-frequency pairings with co- 
existence of the singlet and the triplet states in the case of 
inhomogeneity. 

We consider the general form of the irreducible interaction 
in the path integral (3), 



So = JV„(l)[0 T + f]lMl), 



, 7S (l-2)pJl,2)p yd (l,2), (A-l) 



where the summation is understood for the repeated spin in- 
dices (the greek letters). Due to the hermiticity, the pairing 
interaction satisfies 



Vafiiydi} _ 2) - V 7 g- a p(l - 2)*. 



(A-2) 



Here, the pair field and its complex conjugate have been de- 
fined by 

p a p(l,2) = tf/ a (l)ifrp(2), 

p^(l, 2) = ^,(2)^,(1). (A-3) 

The anti-commutation relation (Pauli principle) is expressed 
as 

p afs {\,2) = - P/te (2, 1), p ap {\,2) = -p pa {2, 1). (A-4) 
By the Stratonovich-Hubbard transformation, we obtain 

S aux = - Yp J j 2 [Pafii 1 > 2 )V( 1 . 2) + V( 1 , 2)p aP ( 1,2)], 
5A = -2^IIf V(1 - 2rl U^ (1 ' 2) ^ (1 ' 2) ' 

(A-5) 

where V(l -2) _1 is the matrix inverse of V{\ -2) whose (a/3)- 
(y6) component is given by V a p- y s(l - 2). Note that the com- 
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plex c-number auxiliary fields A^l, 2) and A^l, 2) satisfy F+(l,2) = -- f f G^ 0) '(l,3)A t (3,4)G fi/? (4,2), 
the Pauli principle as similar to (A 4). PJ3J4 ay yS 

As was explained in the main text, the correct minimum 

(A- 14) 

of the free energy within the saddle-point approximation is where we have introduce d the bare green functions, 

Gf {\,2Y x = -j3 2 <5(1 - 2)5 ap {d T + £), 



obtained by the gap functions that satisfy the gap equations, 
K a p{\,2) = -V aP ; 7S {\-2)F 7S {h2), 

V(l,2) = -V a p- yS (l - 2)^(1,2), (A-6) 



G®hl,2y l = -$ l 6{\ - 2)S aP {d T -&. (A-15) 



Now, let us introduce the fourier transformation according 

and the minimum-energy condition, t0 ^ first ^ of (6)> then wg obtain 

A a p(l,2T = V(l,2) , Al(2, 1). (A-7) ^ ^ = _ (^(^(^ , 

Here we have introduced the anomalous green functions de- ,_ , 

fined by ^ k>) = - (^-k^-k'))^ , 

2) = (^(1)^(2)) mf , F af! (k, k') = (if fa (k)M-k')) MP . 

F a/3 (1, 2) = (^(2)^(1)) mf ee f£,(2, 1), (A-8) *0 = (^(-^(^)) MF = fc). (A-16) 

and A 1 and F + are used for notational simplicity. From the Here, A(£, £') = p~ l jj f 2 A(l,2)e~ ,( - kj:, ~ k ' X2 \ The bare green 

minimum-energy condition (A-7) with the help of the gap functions are given by 

equations, we obtain the relation, „ ( 0)., ,. , r ij , woy, ,/ ,n / A n\ 

_ fe ) = ( l0) n - ft) fe'^QiS = -Gpa (~ k ' ( A ' 17 ) 

F„«(l,2)* = F aB {\,2) = F+ (2,1). (A-9) ^ . , , , , . 

p ' p f ay For the gap functions, we define the founer transformation as 

Integrating out the fermion fields by noting the quadratic ^ , i(kx 1 -k , x 2 ) 

form of the action, ' A a p{\,2) = 2_ l A aP {k,k)e 1 2 , 

kk' 

s +s am = -L l£(paW> WD)^(U)|p| V(i,2) = ^yw) e +'-^) 

M a ,(l,2)-^ M t /12) ^ (1 _ 2)M5t _^J, V 1 ' 2 )-^^*^ ■ < A18 > 

(A 10) r pj ie mm j mum . ener gy conc iition (A-7) and (A-9) become 

we obtain the free energy as ^ ky = ^ ^ = ^ ^ 

Tmf = ~ ff[ V Q ~ 2 r%, yS V(1.2)M1,2) F<*(*,*T = = F^tf,*). (A-19) 

In the case of a homogenous state, we can set F a p(k, k') = 

1 In ( dgt ^ ) . (A-ll) F a p(k)6kk>, A af) (k,k') = A a p(k)6 k k' and so on. Using the 

2[5 \detM / momentum-frequency basis, the free energy (A- 11) becomes 

where M = M(A = A 1 " = 0). 1 

Let us introduce the green functions together with (A-8), Tmf — - ^ ^ Waj8;?tf(& - k')A a p(k)* A y s(k') 

G^(1,2) = -(^(1)^(2)) mf , 



20) 



_ IVi / detM ^M rA 

G^(l,2) = -(^(l)^(2)) MF = -G j3tt (2,l). (A-12) 2/3 ^ " ldetM (*)/ ' 

When the action is quadratic, the matrix M and the green where [V(ry l ] af) - yS = £ 9 W a p;ys{q)e iqr . Using the gap equa- 

functions are related as tion, 

M a p(h2) = ^:^) -Gyhl)} 1 ' (A ' 13) K P {k) = - l -^V al3 , ye {k-k>)F y5 {k'), (A-21) 

and it gives the Gor'kov equations, the free energy is rewritten as 

G*(l,2) = Gg(l,2) - I f rG^(l,3)A y5 (3,4)^(4,2), ^ = J_ V L,^F a , W - lnf^|^-)l . (A-22) 

VJiJ* l P k^ \detM (fe)/_ 

G^g(l,2) = G^'(l,2) - X - f rG^'(l,3)A^(3,4)F tfj8 (4,2), When we decompose the gap function into the singlet com- 

" 3 4 ponent do(k) and the triplet component d(k) as 

F a p(l,2) = -jj j^j^ y (l,3)A yS (3A)G' s/} (4,2), A afs (k) = d Q {k)(icr?) a p + d{k) ■ (/o-ct'V (A-23) 

the explicit expressions are obtained as follows: 

detM(fc) = [col + E + (k) 2 ][to 2 n + E_(k) 2 ], (A-24) 
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with 



D (k) = \d Q (k)\ 2 + \d(k)\ 2 , 

D(k) = d (k)d(ky + d(k)d (ky + i [d(k) x d(k)*] , 



E ± (k) = + D (k) ± \D(k)\. 



(A-25) 



Note that D (k) and D(k) are real functions, and the Pauli 
principle demands, 



do(k) = do(-k) = ±d Q (-k, ico„), 
d(k) = -d(-k) = +d(-k, ia> n ), 



(A-26) 



where the upper and the lower signs correspond to the even- 
frequency and the odd-frequency pairings, respectively. In the 
case of the unitary pairing, i.e., D (k) = 0, no time-reversal 
symmetry breaking occurs even for the odd-frequency pairing 
(it breaks only the /mag/nary-time-reversal symmetry of the 
relative coordinate). 

The green functions are written as 

Gafi(k) = -G'pjrk) = G {k)5 a p + G(k) ■ <r ap , 

F als {k) = = F {k){io- y ) a p + F(k) ■ {Ura^ap, 

(A-27) 

where the scalar and the vector components are given by 



G (k) 
G(k) = 

F o(k) = 

F{k) = 



(ifc> H +ft)[^ + ff + Z>o(fc)] 
[a> 2 „ + E + (k) 2 ][col+E-(k) 2 Y 

(iu) n + £k)D(k) 
[<o 2 n + E + (k) 2 ][oj 2 n + E4k) 2 y 

{oi 2 n + &do(k) + [d Q (k) 2 - d(k) 2 ]d (ky 



[04 + E + (k) 2 ][o4 + E-(k) 2 ] 

(oj 2 „ + f k )d(k) - [d Q (k) 2 - d(k) 2 ]d(ky 



[(4 + E + {ky][u? n + EAk) 2 ] ' (A ' 28) 

For the magnetic-field response, we again introduce the 
minimal coupling in £ of (A- 10), and expand the free energy 
up to 0(A 2 ) by using (A- 13). Then, we obtain 

1 2 "J* 

Tmf(A) =fMF + r^VV K^A^Aji-q), (A-29) 

'j q 



where the Meissner kernel has been introduced by 

K,M) = ^ { G ^)Gp a (k -q) + F ap {k)F+ a {k - q)) 



+ 6ijG aa {k) 



, (A-30) 



and q = (q, ie m ) with the bosonic Matsubara frequency, e m 
2nrnip. 
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